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Abstract 

We report the observation of a near-threshold enhancement in the coJ/ijj invariant mass distri- 
bution for exclusive B — > KujJ/il^ decays. The results are obtained from a 253 fb~^ data sample 
that contains 275 million BB pairs that were collected near the T(4S) resonance with the Belle 
detector at the KEKB asymmetric energy e+e~ collider. The statistical significance of the cjJ/^ 
mass enhancement is estimated to be greater than 8a. 

PACS numbers: 14.40.Gx, 12.39.Mk, 13.25.Hw 
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Recently there has been a revival of interest in the possible existence of mesons with a 

more complex structure than the simple qq bound states of the original quark model. There 

are long-standing predictions of four-quark qqqq meson-meson resonance states ^1] and for 

n 

qq-gluon hybrid states 2\. Searches for these types of particles in systems that include a 
charmed-anticharmed quark pair (cc) are particularly effective because for at least some of 
these cases, the states are expected to have clean experimental signatures as well as relatively 
narrow widths, thereby reducing the possibility of overlap with standard cc mesons. 

B meson decays are a prolific source of cc pairs and the large B meson samples produced 
at 5-factories are providing opportunities to search for missing cc charmonium mesons as 
well as more complex states. From studies of RgK^n^ systems produced in exclusive 
B KKgK~7i~^ decays with a 45 million BB event sample, the Belle group made the 
first observation of the rjc{2S) P]. This state was subsequently seen in two-photon reactions 
by other experiments and in the exclusive e"'"e~ J/iprjc{2S) production process by 
Belle jbl. With a larger sample of 152 million BB events. Belle discovered the X(3872) as 
a narrow peak in the Tr~^ tt~ J/ ip mass spectrum from exclusive B KiT^n~ J/ip decays 0|. 
This observation has been confirmed by other experiments [7]. The properties of the X(3872) 
do not match well to any cc charmonium state 8]. This, together with the close proximity 
of the X(3872) mass with the m^io + m£)*o mass threshold, have led a number of authors to 
interpret the X(3872) as a D^D*^ resonant state |9|. 

In this Letter we report on a study of the uJ/ip system produced in exclusive B — > KuJ/ijj 
decays. We observe a large enhancement in the uJ/ip mass distribution near the uJ/ip mass 
threshold. These results are based on a 253 fb~^ data sample that contains 275 million BB 
pairs collected with the Belle detector. 

The Belle experiment observes B mesons produced by the KEKB asymmetric energy 
e"'"e~ collider [lO|. KEKB operates at the T(4S') resonance {^/s = 10.58 GeV) with a peak 
luminosity of 1.39 x 10^^ cm~^s~^. 

The Belle detector is a large-solid-angle magnetic spectrometer that consists of a silicon 
vertex detector, a 50-layer cylindrical drift chamber (CDC), an array of aerogel threshold 
Cherenkov counters (ACC), a barrel-like arrangement of time-of-flight scintillation counters 
(TOF), and an electromagnetic calorimeter (ECL) comprised of CsI(Tl) crystals located 
inside a superconducting solenoid coil that provides a 1.5 T magnetic field. An iron flux- 
return located outside of the coil is instrumented to detect Kl mesons and to identify muons 
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(KLM). The Belle detector is described in ref. 

We select events of the type B Kn^-n-^n^J/tp, where we use both charged and neutral 
kaons We use the same charged kaon, pion and J/ip requirements as used in ref. 0|. 
For neutral kaons we use tt+tt^ pairs with invariant mass within 15 MeV of rnxg and a 
displaced vertex that is consistent with Kg — > tt+tt^ decay. We identify a vr*^ as a 77 pair 
that fits the 7r° — > 77 hypothesis with < 6. We further require the energy asymmetry 
— i?^2 1/1-^71 + -^72! < 0-9 laboratory-frame momentum to be greater than 

180 MeV. Events with a tt^tt^J/iIj invariant mass within 3a of m^/ are rejected in order to 
eliminate B KTr^ip'; ip' — ^ -K^-K^J/ip decays. The level of e^e~ ^ qq {q = u,d,s or c- 
quark) continuum events in the sample is reduced by the requirements R2 < 0.4, where R2 
is the normalized Fox- Wolfram moment and | cos^^^l < 0.8, where 6b is the polar angle 
of the i?-meson direction in the center-of-mass system (cms). 

At the T(4S'), BB meson pairs are produced with no accompanying particles. As a result, 
each B meson has a total cms energy that is equal to -Ebeam? the cms beam energy. We 



identify B mesons using the beam-constrained mass Mbc = y -^bcam ~ the energy 

difference AE = -Ebcam — Eb, where pb is the vector sum of the cms momenta of the 
B meson decay products and Eb is their cms energy sum. For the final state used in 
this analysis, the experimental resolutions for Mbc and AE are approximately 3 MeV and 
13 MeV, respectively. 

We select events with Mbc > 5.20 GeV and \AE\ < 0.2 GeV for further analysis. For 
events with multiple vr'' entries in this region, we select the 77 combination with the best 

value for the 7r° 77 hypothesis. Multiple entries caused by multiple charged particle 
assignments are small (~ 4%) and are tolerated. The signal region is defined as 5.2725 GeV< 
Mbc <5.2875 GeV and \ AE\ < 0.030 GeV, which is about ±2.5cr from the central values. We 
identify three-pion combinations with 0.760 GeV < M(7r+7r^7r'^) < 0.805 GeV as candidate 
u mesons. To suppress events of the type B — > KxJ/'ip', Kx — > Koj, where Kx denotes 
strange meson resonances such as i^'i(1270), fri(1400), and -ft'|(1430) that are known to 
decay to Kuj, we restrict our analysis to events in the region M{Kuj) > 1.6 GeV. 

Figure ^^a) shows the Mbc distribution for selected events that are in the AE and lj 
signal regions. The curve in the figure is the result of a binned likelihood fit that uses a 
single Gaussian for the signal and an ARGUS function for the background. The fit gives 
a signal yield of 219 ± 23 events. Figure Q^b) shows the AE distribution for events in the 
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FIG. 1: (a) Mbc distributions for B — > KojJ/ijj candidates in the A£^ and uj t:~^t:~t:^ signal 
regions, (b) distribution for events in the and uj signal regions, (c) M(7r+7r~7r^) distri- 
bution for events in the Mbc and /S.E signal regions. The curves are the results of fits described in 
the text; the arrows indicate the signal region for each quantity. 

Mbc and uj signal regions. Here the curve is the result of a fit that represents the signal 
with a single Gaussian and the background with a first-order polynomial. The signal yield is 
196±21 events. Figure^c) shows the M(7r+7r^7r°) distribution for all events in the Mbc-Ai? 
signal region. The peak is well fitted with a Breit-Wigner with the mass and width of the 
C(j(780), broadened by an experimental resolution of 8 MeV. Here the signal yield is 204 ±20 
events. The reasonable consistency in the signal yield from all three distributions indicates 
that KuoJ/ip is the dominant component oi B ^ Kix^ti^-k^J/iI) decays with M(37r) in the 
uo mass range and M{Kuj) > 1.6 GeV. The arrows in the figures indicate the signal regions 
for the plotted quantity. 

Figure El shows the Dalitz plot of M'^{ujJ/ip) (vertical) vs M'^{ujK) (horizontal) for B — >■ 
KujJ/ip candidates in the signal regions. Here the M{Kuj) > 1.6 GeV requirement has been 
relaxed. The clustering of events near the left side of the plot are attributed to B ^ KxJ/'ip', 
Kx — * Kuj decays. There is an additional clustering of events with low uJ/ip invariant 
masses near the bottom of the Dalitz plot. This clustering is the subject of the analysis 
reported here. 

The fits to the Mbc and AE" distributions of Figs. CJa) and (b) indicate that about half of 
the entries with M{ujK) > 1.6 GeV in the Dalitz plot of Fig.|21are background. To determine 
the level of i? ^ Kuj J /if) signal events, we bin the data into 40 MeV- wide bins of M{ujJ/tp) 
and fit for B meson signals. The histograms in Figs. ini^a)-(l) show the Mbc distributions for 
the twelve lowest M{ujJ /ip) bins for events in the AE and uj signal regions. Here there are 
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FIG. 2: Dalitz-plot distribution for B KtoJ/ip candidate events. The dotted line indicates the 
boundary of the M(Kuj) > 1.6 GeV selection requirement. 

distinct B KuJ/ip signals for low uJ/ip invariant mass bins, especially those covered by 
Figs. Efb) and (c). We establish the B — * KuJ/tp signal level for each M{(jjJ/%1)) bin by 
performing binned one-dimensional fits to the Mbc and AE distributions for events in that 
interval using the same signal and background functions that are used to fit the integrated 
distributions of Figs, ^a) and (b). For these fits, the peak positions, resolution values 
and background shape parameters are all fixed at the values that are determined from the 
fits to the integrated distributions, and the areas of the Mbc and Ai? signal functions are 
constrained to be equal. The curves in Fig. El indicate the results of the Mbc fits. 

The i?-meson signal yields from the fits to the individual bins are plotted vs M{ujJ/tp) in 
Figs.El^a) and (b). An enhancement is evident around M{ujJ/ip) = 3940 MeV. The curve in 
Fig. ^a) is the result of a fit with a threshold function of the form /(M) = Aoq*{M), where 
q*{M) is the momentum of the daughter particles in the ujJ/ip restframe. This functional 
form accurately reproduces the threshold behavior of Monte Carlo simulated B KuJ/ijj 
events that are generated uniformly distributed over phase-space. The fit quality to the 
observed data points is poor {x^/d.o.f. = 115/11), indicating a significant deviation from 
phase-space; the integral of f{M) over the first three bins is 16.8 events, where the data 
total is 55.6 events. 

In Fig. Hfb) we show the results of a fit where we include an S'-wave Breit-Wigner (BW) 
function [15^ to represent the enhancement. The fit, which has 'x^/d.o.f. = 15.6/8 (CL = 
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FIG. 3: Mbc distributions for 5 —> K wJ/'i/; candidates in the A^^ signal region for 40 MeV- wide 
uJ/ip invariant mass intervals. The curves are the results of fits described in the text. 
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FIG. 4: — > KloJ/i/j signal yields vs M(LoJ/tp). The curve in (a) indicates the result of a fit that 
includes only a phase-space-like threshold function. The curve in (b) shows the result of a fit that 
includes an S-wave Breit-Wigner resonance term. 



4.8%), yields a Breit-Wigner signal yield of 58 ± 11 events with mass M — 3943 ±11 MeV 
and width F = 87±22 MeV (statistical errors only) .The statistical significance of the signal, 
determined from ln(>Co / J^max) , where >Cmax and Co are the likelihood values for the 
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best-fit and for zero-signal-yield, respectively, is 8.1a. 

The Kuj invariant mass distribution for Mbc-Aii^ signal region events in the region of the 
M{ujJ/ip) enhancement are distributed uniformly across the available phase space and there 
is no evident Kuj mass structure that might be producing the observed mass enhancement 
by a kinematic reflection. Nevertheless, the possibility that different high-mass Ku partial 
waves might interfere in a way that produces some peaking in the uoJ/ip mass distribution 
cannot be ruled out. 

The M(7r+7r^7r^) distributions for different M{ujJ/ip) mass regions exhibit u tt+tt^tt'' 
signals that track the Mbc-Ai? signal yields. The u signal strength is used to infer that 
(90 ± 18)% of the B Kn+n-n^J/ijj events in the M = 3943 MeV enhancement are 
produced via u —>■ tt+tt^tt'^ decays. 

We study potential systematic errors on the yield, mass and width by repeating the 
fits with different signal parameterizations, background shapes and bin sizes. For example, 
when we change the background function to include terms up to third order in q*, the 
yield increases to 75 ± 10 events, the mass changes to 3948 ± 9 MeV, the width changes to 
r = 100 ± 23 MeV and the fit quality improves: x^/d.o.f. = 10.0/6 (CL=12.4%). However, 
the resulting background shape is very different from that of phase-space. For different bin 
sizes, fitting ranges, M{Kuj) requirements, and signal line shapes we see similar variations. 

For the systematic uncertainties we use the largest deviations from the nominal values 
for the different fits. In the following, we assume that all of the Sir systems are due to — >■ 
71^71^71^ decays and include the possibility of a non-resonant contribution in the systematic 
error. This is the main component of the negative-side systematic error; the change in yield 
for different background shapes contributes a positive-side error of comparable size. The 
effects of possible acceptance variation as a function of M{(jjJ/iI)) on the mass and width 
values are found to be negligibly small. 

To determine a branching fraction, we use the BW fit shown in Fig. IH^b) to establish 
the event yield of the observed enhancement. Monte Carlo simulation is used to estimate 
detection efficiencies of 2.4 ± 0.1% and 0.42 ± 0.02% for B K^ujJj^ and K^uJ/iIj, 
respectively. We find a product branching fraction (here we denote the enhancement as 
r(3940)) 

B{B is:r(3940))i3(r(3940) ^ uJ/i)) = (7.1 ± 1.3 ± 3.1) x 10-^ (1) 
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where the second error is systematic. The latter includes uncertainties in the acceptance, 
and the shape of the function used to parameterize the background and the possibility of a 
non-U component of the vr+vr^vr'^ system added in quadrature. Here we have assumed that 
charged and neutral B mesons are produced in equal numbers at the T(45') and they have 
the same branching fractions to the observed enhancement 

In summary, we have observed a strong near-threshold enhancement in the uJ/ip mass 
spectrum in exclusive B —>■ KujJ/ip decays. The enhancement peaks well above threshold 
and is broad |l7^: if treated as an S'-wave BW resonance, we find a mass of 3943± ll(stat) ± 
13(syst) MeV and a total width F = 87 ± 22(stat) ± 26(syst) MeV. It is expected that a cc 
charmonium meson with this mass would dominantly decay to DD and/or DD*; hadronic 
charmonium transitions should have minuscule branching fractions Ji3]. 

The peak mass of the observed enhancement is very similar to that of a peak observed 
by Belle in the J/ip recoil mass spectrum for inclusive e~^e~ J/ipX events near y/s = 
10.56 GeV This latter peak is also seen to decay to DD*, and a search for it in the 
ujj/ip channel is in progress. In addition, we are examining B KDD* decays for a DD* 
component of the enhancement reported here. 

The properties of the observed enhancement are similar to those of some of the cc-gluon 
hybrid charmonium states that were first predicted in 1978 P] and are expected to be 
produced in B meson decays 3]. It has been shown that a general property of these hybrid 
states is that their decays to meson pairs are forbidden or suppressed, and the 
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where D** refers to 



relevant "open charm" threshold is m/j + m^" — 4285 MeV 
the = (0, 1,2)+ charmed mesons. Thus, a hybrid state with a mass equal to that of 
the peak_we observe would have large branching fractions for decays to J/ip or ip' plus light 
hadrons Q]. Moreover, lattice QCD calculations have indicated that partial widths for such 



24| . However, these calculations 



25|, substantially 



decays can be comparable to the width that we measure 
predict masses for these states that are between 4300 and 4500 MeV 
higher than our measured value. 

Acknowledgments 

We thank the KEKB group for the excellent operation of the accelerator, the KEK 
cryogenics group for the efficient operation of the solenoid, and the KEK computer group 



10 



and the Nil for valuable computing and Super-SINET network support. We acknowledge 
support from MEXT and JSPS (Japan); ARC and DEST (Australia); NSFC (contract 
No. 10175071, China); DST (India); the BK21 program of MOEHRD and the CHEP SRC 
program of KOSEF (Korea); KBN (contract No. 2P03B 01324, Poland); MIST (Russia); 
MESS (Slovenia); Swiss NSF; NSC and MOE (Taiwan); and DOE (USA). 



* on leave from Nova Gorica Polytechnic, Nova Gorica 
[1] See, for example, M.B. Voloshin and L.B. Okun, JETP Lett. 23, 333 (1976); M. Bander, 

G.L. Shaw and P. Thomas, Phys. Rev. Lett. 36, 695 (1977); A. De Rujula, H. Georgi and 

S.L. Glashow, Phys. Rev. Lett. 38, 317 (1977); N.A. Tornqvist, Z. Phys. C 61, 525 (1994); 

and A.V. Manohar and M.B. Wise, Nucl. Phys. B 339, 17 (1993). 
[2] D. Horn and J. Mandula, Phys. Rev. D 17, 898 (1978). 
[3] S.K. Choi et al. (Belle Collaboration), Phys. Rev. Lett. 89, 102001 (2002). 
[4] D.M. Asner et al. (CLEO Collaboration), Phys. Rev. Lett. 92, 142001 (2004) and B. Aubert et 

al. (BaBar Collaboration), Phys. Rev. Lett. 92, 142002 (2004). 
[5] K. Abe et al. (Belle Collaboration), Phys. Rev. Lett. 89, 142001 (2002). 
[6] S.K. Choi et al. (Belle Collaboration), Phys. Rev. Lett. 91, 262001 (2003). 
[7] D. Acosta et al. (CDF-II Collaboration), Phys. Rev. Lett. 93, 072001 (2004); V.M. Abazov et 

al, (DO Cohabor ation) , |hep-ex/0405004| Phys. Rev. Lett. 93, 162002 (2004); and B. Aubert et 

al, (BaBar Collaboration), |hep-e x/0406022] 
[8] S.L. Olsen, piep^ex/0407033 and K. Abe et al (Belle Cohaboration), hep-ex/0408TT6l 
[9] N.A. Tornqvist, Phys. Lett. B 590, 209 (2004); E.S. Swanson, Phys. Lett. B 588, 189 (2004); 

F.E. Close and P.R. Page, Phys. Lett. B 578, 119 (2003); S. Pakvasa and M. Suzuki, Phys. 

Lett. B 579, 67 (2004); C.-Y. Wong, Phys. Rev. C 69, 055202 (2004); and E. Braaten and 

M. Kusunoki, Phys. Rev. D 69, 114012 (2004). 
[10] S. Kurokawa and E. Kikutani, Nucl. Instr. and. Meth. A 499, 1 (2003). 

[11] A. Abashian et al (Belle Collaboration), Nucl. Instr. and Meth. A 479, 117 (2002) and 

Y. Ushiroda (Belle SVD2 Group), Nucl. Instr. and Meth. A 511, 6 (2003). 
[12] The inclusion of charge-conjugate modes is always implied. 
[13] G.C. Fox and S. Wolfram, Phys. Rev. Lett. 41, 1581 (1978). 



11 



[14] H. Albrecht et al. (ARGUS Collaboration), Phys. Lett. B 299, 304 (1989). 

[15] We use the form (jv/2_jv/p^[i^((l°r(g»/gp))2 , where Mq is the peak mass and = q*{M = Mq). 

[16] The signal yields for the separate K'^ajJ/tp and Ksujj/tp channels are consistent. 

[17] The Ad{(jjJ mass resolution is ~6 MeV and much narrower than the observed enhancement. 

[18] See, for example, T. Barnes and S. Godfrey, Phys. Rev. D 69, 054008 (2004) and E. Eichten, 

K. Lane and C. Quigg, Phys. Rev. D 69, 094019 (2004). 
[19] P. Pakhlov, |hep-ex/0412042j 

[20] F.E. Close, I. Dunietz, P.R. Page, S. Veseh and H. Yamamoto, Phys. Rev. D 57, 5653 (1987); 

G. Chiladze, A.F. Falk and A.A. Petrov Phys. Rev. D 58, 034013 (1988); and F.E. Close and 

S. Godfrey, Phys. Lett. B 574, 210 (2003). 
[21] N. Isgur, R. Kokoski and J. Baton Phys. Rev. Lett. 54, 869 (1985). 

[22] F.E. Close and P.R. Page, Nucl. Phys. B 443, 233 (1995); F.E. Close and P.R. Page, Phys. 

Lett. B 366, 323 (1996). 
[23] F.E. Close, Phys. Lett. B 342, 369 (1995). 

[24] C. McNeile, C. Michael and P. Pennanen (UKQCD Collaboration), Phys. Rev. D 65, 094505 
(2002). 

[25] C. Banner et al, Phys. Rev. D 56, 7039 (1997); Z.-H. Mei and X.-Q. Luo, Int. J. Mod. 
Phys. A 18, 5713 (2003); and X. Liao and T. Manke, plep^lat/0210030| 



12 



